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Abstract 

We use multiscale-multispace correlations and Fourier transform techniques, to study some 
intermittent random field properties, which escape analysis by structure function scaling. These 
properties are parametrized in terms of a set of scale ratios, giving the typical interaction distances 
in space and scale of the random field fluctuations, and the characteristic lengths over which 
these fluctuations act coherently to generate intermittency. The relevance of these techniques in 
turbulence theory is discussed. 
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Turbulence is usually depicted as a situation of irregular fluid motion, with energy being trans- 
ferred from large to small scales, by interaction between eddies of increasingly smaller size. This 
mechanism of energy transfer, usually referred to with the name of cascade, takes different forms 
depending on the choice of turbulence description. In closure Q| , a Fourier representation of the 
Navier-Stokes equation is typically used. In this case, the transfer of energy between eddies is 
expressed in terms of an energy flux in wavevector space. This flux is a nonlinear kernel of the 
spectral energy density, and therefore a fully global quantity, carrying no information on spatial 
structure and transfer fluctuations. In intermittency modelling, instead, the description is neces- 
sarily local in space, at least implicitly. In the random beta model ||, for example, this is obtained 
using a superposition of nested eddies of increasingly smaller size, with a fluctuating eddy to eddy 
energy transfer. 

In both circumstances, the dynamics forces correlations between the fluctuations at different 
scale and space locations. Only recently, however, there has been some interest in studying cor- 
relations in turbulence ||, [Q, beyond the standard approach of focusing on structure function 
scaling. 

A tool that is receiving great attention, especially in the study of intermittency, is synthetic 
turbulence ||, @, @- Synthetic turbulence has no pretense of modelling turbulence dynamics, its 
goal being more the one of reproducing the statistics observed in experimental turbulence data. 
Even so, the multiplicative noise algorithms used to produce the turbulent signal contain implicit 
assumptions on the nature of the turbulent flow that go beyond the task of reproducing anomalous 
structure function scaling. 

Typically, these synthetic signals are expressed as a superposition of wavelets with random 
coefficients: 

N 

^(x) ^^2Y, A ^ W s( k h n ,Vh n ,x), (1) 
n h n 

where ^(x) is the synthetic turbulent signal, with ws(k,y,x), wavelets of scale k and position y 
||, and the vector index h„ = (ho, hi, ....h n ) identifying the position of the wavelet in the cascade 
through the sequence of its ancestors: the integer h n labels the h n -th daughter wavelet generated 
at the n-th step in the cascade, by wavelet h n _i. The standard practice || is to arrange the 
wavelets on a tree structure in k — y space, with fch„ = k n = k~o2 n , hi — ±1 yh n — ?/h„_i + ^f-- 
The amplitudes Ah n are then generated through a random multiplicative process. 

This picture closely parallels the one of a random beta model the resulting correlation 
pattern has clearly a multifractal (multiaffine) nature, depending on the singularity (lack of singu- 
larity) of ^(x), as N — > oo in Eqn. (1). The difference, in the case of the beta model, is that the 
multiplicative random process, arises explicitly out of the assumption of an energy transfer, local 
in both space and scale. 

The purpose of this letter is to present some techniques, based on the use of wavelet and Fourier 
analysis, to extract informations on the algorithm of signal generation, directly from the signal 
statistics. We focus in particular on two important aspects: the determination of which "building 
block" wavelet ws is being used in the synthetic signal, and which is the typical interaction distance 
between wavelets, as expressed by the separation in k — y space: (fch„ — feh„_i,2/h„ — 2/h„_i)- In 
real turbulence, these aspects have a direct interpretation in terms of eddy structure and energy 
transfer properties. 

The choice of the wavelet, identifies a cell of k — y space, which acts coherently to generate 
intermittency; this is in practice determined by the "number of wiggles" of the wavelet, i.e. the 
product as = kXs of its characteristic wavevector and spatial extension. We can thus focus on the 
case in which ws is a Gaussian wavepacket, without fear of loosing to much in generality: 

ws(k, y, x) = exp(ik(x - y) - (x - y) 2 / A|). (2) 

The interaction distance provides informations on the degree of non-locality of the dynamics, by 
telling us whether only nearest neighbor cells in k — y space can interact directly, or longer range 
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interactions are allowed. The algorithms of turbulence synthesis studied in ||, || and j| are strictly 
local. To allow for non-locality, the simplest approach is to treat the wavelet coordinates in k — y 
space, on the same ground as the amplitudes A^ n , as stochastic variables. In this way, the wavelets 
take more the aspect of eddies, with complete freedom of location (and overlap) in k — y space ||. 
The dynamics is described through the probability for the transition (rn^4h„_i,hifch„_i,?/h„_i) ~> 
(In Ah n , In kh n , j/h„), which is taken in the form: pa(A' /A\k' /k)pk(k' /k)p x (k\y— y'\), with J d\nyy p 
xpA(y\x) = c p x~^ p to insure power law scaling in k for the moments of A. The phases are supposed 
random so that intermittency in ^ arises simply as a consequence of the intermittency of the 
amplitudes A. We take for simplicity, p x and pk to be Gaussians respectively in space separation 
and In /ch„/^h„_i ■ The transition probabilities over n cascade steps will be Gaussians as well: 

k ( k 2 \y-y'* 2 ^ 



P x [n, k\y -y\)zt —rj — exp I j— J 
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P k (n; m = -4- exp ( - Mk'/k)- n^ 

with b giving the amount of space non- locality and Az/z the degree of discreteness in scale, in 
the cascade algorithm. Thus, the wavelet wavevectors in the cascade will be centered at k n = 
koexp(nz) with increasing width Ak n ~ n^Azko- The probability pa, which is the multiplier 
distribution for the coefhncients Ah n , can be obtained starting from the distribution of the scaling 
exponents Q q , using standard techniques p0| . We do not need its explicit form here, however. 

To measure effects such as the degrees of interaction non-locality and discreteness, in an in- 
termittent signal of unknown origin, multiscale-multispace wavelet correlations would be clearly 
the instrument of choice Q. This is conditioned, however, to having some knowledge of the 
structure of the building block wavelets ws, especially, as regards the value of the parameter as- 

The simplest wavelet correlations are (IvE'fcj/l 2 ) and ky\ 2 \^ k 1 y'\ 2 ) , where ^ky is the component 
of W on a wavelet (—k)~ s d x WA(k, y, x). Here, wa is taken again to be a Gaussian wavepacket, with 
a a — k\A not necessarily equal to as, and s is chosen large enough to kill large scale contributions 
to *: 

H> ky = X A x k- s [ dxw A (k,y,x)d s x y(x). (4) 



In a cascade model like the one considered here, the wavelet components ^>ky entering (l^tj,! 2 
x Y&k'y' | 2 ) ma y come from a single building block wavelet ws ("one-eddy" contribution), if |fc — k'\ 
and \y — y'\ are both small enough; otherwise they will come predominantly from two distinct ws 
placed at (k,y) and {k',y') ("two-eddy" contribution). If \y — y'\ < L, with L the largest scales 
in the signal, there will be a correlation between the two building block wavelets ws, due to their 
coming from some common ancestor at scale k ~ bas\y — This situation is expressed through 
the formula: 

(\A hn \ 2 \A h ,J 2 ) =c i (k h , m L)-^(k h Jk h ,J-^(k hL /k hp )^- 2C * (5) 

with fch„ > kh^ and p the cascade step at which the genealogical tree of h„ and h' m branches: 
hi = h! i for < % < p and hi ^ h^ for i > p. Thus, the lower the branching takes place in the tree, 
the closer the correlation gets to its disconnected limit (\Ah n | 2 |-Ah^ | 2 ) ~ (kh' m kh n )~ ( ' 2 ■ 

Indicate with C(k, y; k', y') the square modulus of the component of ws with respect to wa- 

C(k A ,yA;k s ,ys) 



■ exp 



A a 1 fc A S / dxw A( k A,yA,x)d x w s (k s ,ys,x) 
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where Ak = k A — k$ and Ay — y A — ys- We can now write down the expression for the second 
order correlation. From Eqn. (3), for InfcL > z 3 /Az 2 , discreteness effects can be neglected, and we 
have, indicating with (Pk(n,k/ko)), average over the size of the large scale building block wavelets: 

(\*ky\ 2 )=c 2 f dydlnfcVfc„(P fc (n,fc/fco))(fci)- C2 C(fc,y;fc,y)~^^(fcL)-^. (7) 
J ^ za-A 

n 

where the difference e = ( 2 — (2 = 0{z~ 3 Az 2 ) is due to the fluctuation in the scale ratio between a 
daughter wavelet h„ and her parent h„_i, and we have used the result, valid for lnfc/fc > z 3 /Az 2 : 
^^ 1 (fc„/A:)Pfc(n, k/ko) — z^ 1 (k/k ) e . We have a similar expression for the one-eddy contribution 
to the 4-th order correlation: 

<l**„| 2 |*Wl 2 >i ^ ^^(kL)-^(k'/k)-^C(k,y;k',y>) (8) 
za A 

From Eqns. (8) and (6), if the space-scale separation is large enough, (l^fei/Pl^fe'y'l 2 } wm be 
dominated by the two-eddy contribution. In this case, however, it is not sufficient to impose that 
kL be large, to avoid discreteness effects. Limiting the calculation to the connected part of the 
correlation, we find: 



p pin k 

\^ky\ 2 \^k'y'\ 2 )c = 2c 4 / dlnfcdlnfc'dydy' / dlnfc {kL)-^(k' /k)-^(k/k) u - 2C - 2 

Jk'>k JliiL- 1 

k n km 



x E E E ^( P *(P> k/h))Pk(n - p, k/k)P k {m - p, k'/k) 

n m— Op— 1 P 

xP x (k\y - y'\)C(k, y- k, y)C(k', y'; k', y'). (9) 

where (k m /k p )P x (k\y — y'\) is the space density at y', of wavelets h' m generated from the branching 
at hp, given the presence of a wavelet h„ at y, and k n is the space density of wavelets h„. From here, 
informations on the space-scale structure of the intermittent signal can in principle be obtained. 
The slow decay of correlations between wavelets, shown in Eqn. (5), however, suggests that the 
effect of coherency, important for k\y — y'\ < as, and the one of interactions, important for 
k\y — y'\ < bas, will be superimposed in a way that is difficult to disentangle. This warns against 
trying to measure as looking for crossovers in the space scaling of correlations. We have in fact 
two mechanisms of correlation distruction: one is the decay of the one-eddy contribution, due to 
the two analyzing wavelets ceasing to overlap with a single u>s- The second is the decay of the 
two-eddy contribution, i.e. the decrease with distance of the correlation between different u>s- 
Only the first effect, of course, gives informations about the shape of ws- 

Surprisingly, this strongly local information is more easily obtained using a fully global in- 
strument like a Fourier transform, at least if the random phase hypothesis is satisfied. The rea- 
son is the maximal definition in scale of a Fourier component. In this case, a correlation like 
(\E r2 _ fe , 3/fe— A^fc+zx), with standard Fourier components, will be non-zero only when A < A^ 1 = 
k/as', in fact, a 4-th order correlation can receive contributions at most by two ws (plus their 
complex conjugate), which cannot cover a section of k including at the same time —k, k — A and 
k + A, when A is large. In the example of our cascade model, we can calculate this fourth order 
correlation rather easily. In the case of negligible discreteness effects, we find: 

where a = (2 + ( 2 + a |/2)^ and the two terms in square brackets are respectively the one- and 
two-eddy contribution to the correlation. 
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Once as is known, it is possible to set a a = as in the analyzing wavelet, in order to maximize the 
overlap between the analyzing wavelets wa and the coherent regions identified by the building block 
wavelets Ws- At this point it is possible to look at the dependence on space and scale separation, 
of the expression for the correlation ky\ 2 \^ k'y'\ 2 )c provided by Eqn. (9). We consider first the 
case in which discreteness effects can be neglected, corresponding to the regime Ink' /k > z 3 /Az 2 . 
Setting cla = as we obtain: 



<l^l 2 l* fc vl 2 K 



2 -^l {k L)-Hk'/k) 

z A ba s Jo 



In kL 



x exp 



(C 4 -2^)*-(e 2 * + F) 



"'"L-'/L^ / dx(l + Ae- 2x 
ifc 2 



\y - y 



t\2 



(11) 



where F = b (1 + (k/k') I ). From inspection of this equation, we see that the integral receives 



contribution, for bas large enough, from max(0, In 



k\y-y'\ ■ 



, InF) < x < lnfcL, where the integrand 



is essentially exp ((£4 — 2^)x). We find then, for 2^ — (4 small: 

2C47T^ 



(l*^| 2 |*fe' y '| 2 ) c * -^(kL)-^(k'/k)-^ max (uGkL, (2& - Q^G 2 ^ (12) 



where G 



z 3 bas 

max(l, k \y~ y I , F i). We thus have a close range and a large separation range, with 



transition at \y — y'\ 
121 



(as/k) max(6, (1 

2\ 



{k/k' 



only in k: {\^ky\ 2 \^k' y '\ 2 )c ~ (kL) ^(k'/k) <* 2 lnkL, while, in the large separation regime, the 



In the close range, the correlation scales 

n regim 

'|C4-2C 2 



correlation contains a factor scaling like a power in \y — y'\: (\^ky\ \^k'y'\ )c oc \y — y'\^~' A ^. This 
slow decay in \y — y'\, observed also in ||, and the one of structure functions, with respect to a a 
||, have the same origin in the dependence on fch p of (|Ah„ | 2 |^4h^ | 2 ) [see Eqn. (5)]: larger space 
separations between building block wavelets imply a lower branching in their genealogical tree. 

Thus, if b is large, for any scale k there is a second characteristic length, beyond As, which 
gives the degree of space non-locality in the interaction, and which can be measured by locating 
the crossover to power law scaling with respect to \y — y'\, in (|^fci/| 2 |v& r k'i/'| 2 )- 

If z 3 / Az 2 is large enough, there is a range of scale separations: lnfc'/fc < z 3 / Az 2 where 
discreteness effects are relevant. In this regime, only one term contributes in the sums over m and 
n in Eqn. (9): the one corresponding to the closest k n and k m respectively to k and k' . (The sum 
over p remains continuous, provided kL is large). Setting again a a = as, the result, after some 
lengthy algebra, is equivalent to the one of Eqn. (11): 



I* 



2 l**vl 2 > 



2C47T2 z 2 

Azbas 



In kL 



dxH-^(l + Fe- 2x ^ 



x exp 



(C 4 -2Qx-(e 2 * + F) 



-ifc 2 



y-y 



/|2 



b 2 a 2 



2HAz 



(13) 



where H = 2x + lnfc'/fc 



of 



Az 2 a A s 



and Akk' = z 1 lnfc'/fc — int(z 1 lnfc'/fc) is the decimal part 



lnfc'/fc. The main difference from Eqn. (11) lies in the quadratic term in A^k' in the 



exponential. As in the case of Eqn. (11), the integral receives contribution for InG < x < lnfcL. 
The integrand has a saddle at: x ~ ag(z/2) 2 (2^2 — C4) -3 \&kk> I — 2z 2 which starts playing a role, 
however, only at very small values of Akk 1 ■ Thus the integral can be estimated almost always with 
steepest descent at x ~ InG. The result is: 



I* 



ky\ 



I* 



k'y' 



2C47T 2^2 

Azbas 



(kL)-^(k'/k)-^G 2 ^-^ 



5 



x max (inGfcL, (2fi - Ci)" 1 J exp ( - - j (14) 

which differs from Eqn. (12), again because of the term in A^. This produces oscillations in 
the correlation dependence on lnfc'/fc, (lacunarity), which have period z, and die when Ink'/k or 
In k \ y ^~ y become larger than z 3 Az~ 2 . This at least in the physically interesting regime where the 
uncertainty in the wavevector over a cascade step is of the same order or larger than the spectral 
width of ws- 

Dynamically, what happens is that the correlation between ws at a scale separation In k' jk 
which is not an integer number of z, is due to a common ancestor at a scale k, which must be 
distant a sufficient number n of cascade steps from k. This is necessary for the width n'Az to 
cover the difference Akk' between z~ 1 Ink 1 /k and the next integer. Clearly, the greater the number 
of cascade steps, the smaller the contribution to the correlation, and this explains the decrease in 
A fcfc ,, of (|1'fc2,| 2 |1'fcv| 2 )c, shown in Eqns. (13-14). 

We thus have three additional characteristic parameters describing the structure of a cascade 
generated random field: the oscillation period in Ink' /k of the correlation, and the separation in 
space and scale over which these oscillations die off. As regards turbulence modelling, it is worth 
mentioning, that different values of z correspond to different degrees of non-locality in the Fourier 
structure of the Navier Stokes equation energy transfer terms. Analysis carried on by means of 
direct numerical simulations and closure has shown indeed that such non-local effects may be an 



important component of turbulence dynamics 1 1 1 



Summarizing, for any fluctuation scale k we have up to five additional scales at our disposal, to 
characterize an intermittent random field. This beyond the standard spectrum of scaling exponents, 
provided by structure function analysis. It is worth considering that these quantities have been 
introduced through an operative definition, which allows their measurement independently of the 
signal being synthetic, and originating from a random cascade. In particular, even in the case 
of a generic random field, in which one would expect the coexistence of fluctuations of different 
shapes, a concept like the coherence lenght As is going to maintain its physical meaning, at least 
on an average sense. The same holds for all the quantities obtained from the multiscale-multispace 
correlation (|^ , fej / | 2 |^ , fe' y ' | 2 ): the degree of space non-locality, identified by the crossover to scaling 
in the space separation, the lacunarity period, measured through the oscillation in Ink' /k, and the 
separation in scale and space for the decay of these oscillations. This suggests, that methods for 
the characterization of intermittent signals, like the ones presented here, could be of some interest 
also for people with access to experimental high Reynolds numbers turbulence data. 

From the point of view of turbulence modelling, it is interesting to examine in more detail what 
happens in some limiting cases. 

A first limit is obtained when b — L/\$, and corresponds to the maximum degree of space non- 
locality in the interaction. In this regime, F ~ (As/L) 2 and {\^ky\ 2 \^k'y'\ 2 )c ~ (fcfc') ^ which 
means that the wavelets ws are essentially uncorrelated, although intermittently distributed; in 
fact, the one-eddy contribution leads still to the anomalous structure function scaling: (l^fcj 4 ) ~ 
(fcL) -1 ' 4 . This situation corresponds to the random eddy model, studied in JO). 

A different, more interesting kind of space non-locality is obtained when As ~ L, which implies 
that as = k\s — > oo as k — > oo. In this case, the building blocks of the random field are 
not wavelets anymore, but Fourier components. We would have then a signal generated by a 
superposition of random phase Fourier modes, analogous to the turbulence picture one gets out 
of statistical closures. The only difference would be the increasing intermittency of the Fourier 
amplitudes as k — > oo. However, as we could expect, this does not produce any intermittency in the 
signal; in fact, it can be shown, by generalizing the result of Eqns. (8), (12) and (14), to the case 
of a,A fixed and as large, that: (\^ k y \ 2 )~ 2 (l^ k y \ 4 ) c K a^ 1 - Hence, the non-Gaussian contribution 
to the correlation decays like k~ x as k — > oo. 
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